NONEXTENSIVE THERMOSTATISTICAL APPROACH TO 
THE THERMOLUMINESCENCE DECAY 

0. Kayacan a '* and N. Can a 
a Department of Physics, Faculty of Art and Science, Celal Bayar 

University, 45030 
Muradiye, Manisa-TURKEY 

Abstract 

In this study, thermoluminescence decay is investigated within Tsallis ther- 
mostatistics(TT). We belive that this is the first attempt to handle themolu- 
minescence decay process within TT. 
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I. INTRODUCTION 



The luminescence intensity of a pre-irradiated phosphor is recorded as a function of the rising 
temperature. The resulting curve between the luminescence intensity and temperature, 
known as the glow curve, first increases and then, after reaching a maximum, decreases. If 
n represents the ocuupied trap density at time t, then the following equation gives the rate 
of thermal release of electron at temperature T: 

(111 

^ = -nsexp(-E/k B T). (1) 

If the phosphor is heated at a constant rate of increase of temperature, i.e. dT/dt = 
cons tan t = r, the above equation can be rewritten as 

— = --exp(-E/k B T). (2) 
n r 

If assuming initially at the start of heating, T = T iy n = n , then the integration of 
Eq.(2) results in 



n 



= n exp |~ j\xp(-E/k B T')'dr 



(3) 



The thermoluminescence intensity / as a function of temperature is given as: 



I = -c ( ^ J = cn s [exp{-E/k B T)] 



exp (~ £ r exp(-E/k B T'ydT' 



(4) 



II. TSALLIS THERMOSTATISTICS 

Boltzmann-Gibbs statistics is used to study the systems having the following conditions: 

(i) the spatial range of the microscopic interactions are short-ranged, 

(ii) the time range of the microscopic memory is short-ranged, 

(iii) the system evolves in a Euclidean-like space-time. 
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These kinds of systems are said to be extensive. If a system does not obey these restric- 
tions, Boltzmann-Gibbs statistics seems to be inappropriate and a non-extensive formalism 
must be used. Tsallis thermostatistics (TT) is one of these formalisms. 

TT has been applied to some concepts of thermostatistics [1] and also achieved in solving 
some physical systems, where BG statistics is known to fail: stellar polytropes [2], Levy-like 
anomalous diffusions [3-7], two-dimensional Euler turbulence [8], solar neutrino problem [9], 
velocity distributions of galaxy clusters [10], fully developed turbulence [11-14], electron- 
positron and other high-energy collisions [15-19], anomalous diffusion of Hydra viridissima 
[20] and nematic liquid crystals [21]. 

Since the paper by Tsallis [22], there has been a growing tendency to nonextensive 
statistical formalism. It has been shown that this formalism is useful, because it provides 
a suitable theoretical tool to explain some of the experimental situations, where standard 
thermostatistics seems to fail, due to the presence of long-range interactions, or long-range 
memory effects, or multi-fractal space-time constraints. 

TT has been applied to various concepts of thermostatistics and achieved in solving some 
physical systems, where Boltzmann-Gibbs statistics is known to fail. Recently, Oliveira 
et al. [23] have pursued the idea, based on novel experimental and theoretical results, 
that manganites are magneticallynon-extensive objects. This property appears in systems 
where long-range interactions and/or fractality exist, and such features have been invoked 
in recent models of manganites, as well as in the interpretation of experimental results. 
Therefore we could expect that this kind of materials could be examined within a non- 
extensive thermostatistics, one of which is called Tsallis thermostatistics (TT) and as you 
can see in our paper, TT has been applied to various concepts of thermostatistics and 
achieved in solving some physical systems, where Boltzmann-Gibbs statistics is known to 
fail. On the other hand, we are still working on application of TT on thermoluminescence 
process, so another studies on this subject are forthcoming, experimentally and theoretically. 

TT considers three possible choices for the form of a nonextensive expectation value. 
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These choices have been studied in [24] and applied to two systems; the classical harmonic 
oscillator and the quantum harmonic oscillator. In that study, Tsallis et al. studied three 
different alternatives for the internal energy constraint. The first choice is the conventional 
one, 

w 



The second choice is given by 



W 



(5) 



(6) 



i=l 



and regarded as the canonical one. Both of these choices have been applied to many different 
systems in the last years [25]. However both of them have undesirable difficulties. The third 
choice for the internal energy constraint is 

T w n q F- 



2^i=i Pi 



u (3) 

U q . 



(7) 



This choice is commonly considered to study physical systems because it is the most appro- 
priate one, and is denoted as the Tsallis-Mendes-Plastino (TMP) choice, q index is called 
the entropic index and comes from the entropy definition, 

1 _ V w r) q 
, 1 2^i=i Pi 



S n — k- 



q-l 



(8) 



where A; is a constant, J2iPi — 1 is the probability of the system in the i microstate, 
W is the total number of configurations. In the limit q — > 1, the entropy reduces to the 
well-known Boltzmann-Gibbs (Shannon) entropy. 
The optimization of S q leads to 



P 



(3) 



i-a-^-c/fVE^^f) 9 
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K3) 



(9) 



with 



w 

zf ] = E 



w 



i-(i-?)0te-f{ 3) )/E(p 



(3)^ 
3 I 
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(10) 



This equation is an implicit one for the probabilities p«. Therefore the normalized q 
expectation value of an observable is defined as 

S T W rfl A 



(ii) 



T w v q 
l^i=i Pi 

where A denotes any observable quantity which commutes with the Hamiltonian. This 
expectation value recovers the conventional expectation one when q — 1. As mentioned 
above, Eq.(ll) is an implicit one and in order to solve this equation, Tsallis et al. suggest 
two different approaches; "iterative procedure" and " (3 — > f3' n transformation. 

The luminescence intensity of a pre-irradiated phosphor is recorded as a function of the 
rising temperature. The resulting curve between the luminescence intensity and tempera- 
ture, known as the glow curve, first increases and then, after reaching a maximum, decreases. 
If n represents the ocuupied trap density at time t, thenthe following equation gives the rate 
of thermal release of electron at temperature T within Tsallis thermostatistics: 

dn 



dt 



-ns{\ - {l-q)E/kT) { Hl 



(12) 



It is important to note that the rate is expressed as a power law rather than exponential 
one. If the phosphor is heated at a constant rate of increase of temperature, i.e. dT/dt = 
constant = r, the above equation can be rewritten as 

— = --{l-{l- q )E/kT)^\ (13) 

71/ T' 

We assume that initially at the start of heating, T = Tj, n = n . Then the integration 
of Eq.(13) results in 



n = no 



1 + (!-?) 



V 



H " - q)E/kT>)^dT>) ^ - 1 



V 



1-q 



/ 
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(14) 



The thermoluminescence intensity / as a function of temperature is given as: 

' dn" 
~dt . 



I = —c 
( 



cn s 



! + (!-<?) 



V 



1 - (1 - q)E/kT) ( — > 
f (- s - - (1 - q)E/kT') ^dT') lk -l \ \ 1 



1-q 



(15) 



Figure shows the computed glow curves corresponding to Eq.(15). The parameters in 
the Eq.(15) are: q = 5i^.s _1 , electron trap depth E = 0.8 eV and s = 5 x 10 11 s _1 . In the 
figure, it is shown the effect of the nonextensivity. In other words, if the rate of thermal 
release of electrons at temperature T is expressed as a power law relation, then the glow 
intensity vs temperature behaves as shown in the figure. An interesting point is that in all 
the curves corresponding to different values of q entropic index, the intensity approximately 
has the same value at T ~ 390i^. 

III. RESULTS AND DISCUSSION 

Presumably there is some normalization at 390 K in order to compare with q = 1 
case. The peak sequence moves smoothly with q value, this is in much the same way 
as caused by changing thermoluminescence kinetics if one varies the relative number of 
trapping and recombination sites and includes secondary, or back reactions. Whilst the 
data presented here are of a preliminary nature, it is sufficient to give encouragement for 
further developments. 
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FIGURE CAPTIONS 

Figure. The computed glow curves corresponding to Eq.(15) for various values of q. 



